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Abstract 

We consider a system of three surfaces, graphs over a bounded do- 
main in R^, intersecting along a time-dependent curve and moving by 
mean curvature while preserving the pairwise angles at the curve of 
intersection (equal to 27r/3.) For the corresponding two-dimensional 
parabolic free boundary problem we prove short-time existence of clas- 
sical solutions (in parabolic Holder spaces), for sufficiently regular ini- 
tial data satisfying a compatibility condition. 



0. Introduction: Triple junctions of graphs. 

The goal of this paper is to prove short-time existence for the following 
geometric evolution problem: a configuration of three surfaces in M'^, bor- 
dered by a common time-dependent curve of intersection A(t), moves with 
normal velocity given by mean curvature (at each interior point) in such a 
way that the pairwise angles defined by their unit normals along A(t) are 
constant throughout the evolution, equal to 27r/3 radians. 

The corresponding evolution for systems of curves was considered in 
[BronsardReiticli] (where short-time existence was proved) and in [Mantegazza et ah] , 
which includes a continuation criterion and the blow-up analysis at the first 
singular time. Very recently, the preprint [LensSeminar] includes global 
existence results for symmetric systems of curves of 'lens type'. In addi- 
tion to being the natural time-dependent version of the classical problem 
of minimal surfaces meeting along a 'liquid edge' ( [Dierkes et al.] p. 299 
ff.), the problem derives its interest from being the sharp-interface limit of 
a well-known class of parabolic evolutionary models in materials science, 
defined by a vector-valued order parameter and a multi-well potential. (See 
[BronsardReiticli] for a discussion in the case of curves.) 

In the case of surfaces, the parametric approach adopted to prove lo- 
cal existence for curves in [BronsardReiticli] and [Mantegazza et ah] does 
not work, that is, does not lead directly to a well-posed parabolic system. 
The essential difficulty can be traced back to the fact that the junction is 
now a one-dimensional object, and one does not expect that independently 
evolving parametrizations of each surface in a triple junction configuration 
will continue to match pointwise along the junction, even if they do so at 
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t = 0. On the other hand, it is difficult to encode the matching condition 
analytically so as to allow for such 'sliding' along the junction. Thus one 
is naturally led to consider the special case of graphs, for the simple reason 
that the graph parametrization is canonical, so that the surfaces' meeting 
along A{t) corresponds to three functions coinciding pointwise along a mov- 
ing interface Tf in M^. 

In this way the problem for graphs becomes a more-or-less classical free 
boundary problem for a quasilinear parabolic system in two dimensions. By 
analogy with the usual Stefan problem terminology, this would be a 'three- 
phase problem', with two of the 'phases' coexisting in the same domain. The 
corresponding 'one-phase' problem- the mean curvature motion of a single 
graph over a time-dependent domain D{t) in M", intersecting M" at a con- 
stant angle- has been dealt with in a recent preprint of the author ( [Freirej ) , 
including local existence and some results on long-time behavior. In [Freirej . 
we introduced additional 'orthogonality' conditions at the free boundary to 
obtain a well-posed system. If one tries to extend this method to the more 
general case considered here, compositions of the unknown functions with 
the (also unknown) diffeomorphisms pulling back the evolution to a fixed do- 
main come into play, and since compositions behave poorly in Holder spaces 
it becomes difficult to implement a fixed-point scheme using this method. 
Instead, we use here an ingenious transformation method which we learned 
from the paper [BaconneauLunardi] . (These authors considered a one-phase 
multidimensional free-boundary problem for a semilinear parabolic system) . 

Now for a precise description of the problem and results. We focus in 
this paper on configurations parametrized by two disks and one annulus, 
but it is clear the method would work for more general configurations (for 
example, two annuli over one sub-domain, one disk over the other; or three 
annuli.) In addition, the system lives in a fixed bounded domain 17 C M^, 
and we must specify how one of the surfaces intersects d^l; we choose or- 
thogonal intersection, which leads to a standard Neumann condition at the 
fixed boundary. More generally, one could prescribe orthogonal intersection 
with a non- vertical, fixed 'support surface'. This would introduce additional 
moving boundaries into the problem but analytically would present no es- 
sential difficulty, so we avoided this in the interest of simplicity. 

A configuration at time t consists of three surfaces: 

= graphiwl), wf = w'{-,t) : D{ ^ R, 1 = 1, 2, 3, 

where D] = and = — D] are bounded, connected open subsets of 
the bounded domain J7 C M^. The intersect along an embedded closed 
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curve At C M^, a graph over the simple closed planar curve Tt = . 
The upper unit normals are: 

= ^[-Dw^, 1], = ^Jl + \Dw^\'^. 

Along Ff, we impose the following conditions: 

(i) (matching) w^{y,t) = w'^{y,t) = w^{y,t), yeVt (5(7°). 

(ii) (angle) N^ + N^ = N^. 

That is, the intersect making pairwise 27r/3 angles. In addition, 
intersects the vertical cylinder over orthogonally: if uj denotes the inner 
unit normal to dil, we require: 

da;W^ = on dn. (BC^s) 

The angle condition is equivalent to two scalar conditions on T^: 

f ^ + ^ = (-SC^) 

where n is the inner unit normal of Fj. 

To specify a configuration, we need three time-dependent functions {t) : 
^R, defining graph parametrizations of the E^: 

-.E^^ R^, = {{x, t) G X [0, T];xe 

G^{x,t) = [x,w^{x,t)] e X M. 

(The labeling of the surfaces is chosen so that near A(t), S2 lies below S^.) 

We wish to evolve a given initial configuration through configurations, 
so that on each surface the normal velocity is given by the mean curvature 
at each point: 

{dtG',N') = H'. 

Recall = tVgih^ = g^^^h^j, where is the second fundamental form of 
S^, pulled back to via G: 

h\e,e) = -{deN^,deG^) = {d^G^ {e,e),N^) = ^d^w\e,e), e G M^, 

and is the induced metric: 

gf. = {d,G', djG') = 5ij + diw'djw'. 
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Mean curvature motion of graphs is given by the equation: 

{MCM)^ Lgi [w'\ = 0, Lgi [/] := dtf - iVgiD^f. 



Our local existence result follows: 

Theorem. Let $7 C 6e a hounded open set with C^+" boundary, 
Dl = Dl and = Q - open subsets of Q with common boundary 
component Tq, a simple closed C^"*"" curve in O. Given G C^+'^{DI), 
assume: 

(1) The Sq = graph{wQ) define a triple junction configuration; 

(2) The mean curvatures Hq of Hq satisfy at points o/Fq the following 
compatibility condition: 

Hi + Hi = Hi ye To. 

Then there exists T > depending only on the initial data, domains 
E^cQx [0,T] and functions G C2+°'i+"/2(£;|,) such that: 

(3) is a classical solution of Lgi[w^] = in ; 

(4) The conjugation conditions {BCyj) hold on = dDl, and the bound- 
ary condition (BCyjs) holds on dil. 

Notation. 

E^= [j Di X {t} C n X [0,T]; 
te[o,T] 

it is always assumed that T = sup{t G [0,T);Dl 7^ 0V/}. We use standard 
notation for parabolic Holder spaces; for example, the norm in C'^'*"°'^+°/^(£') 
(where E cQx [0,T]) is given by: 

with the standard notations for Holder-type difference quotients in x or in 
t. a is an arbitrary constant in the open interval (0, 1), fixed throughout 
the paper. 

The subscript t always denotes a map (or function, or set) at time t, 
never a partial derivative (denoted dt)- The unit vectors tq are the inner 
unit normal (resp. counterclockwise unit tangent) vectors to Fq, while n 
and r denote the corresponding unit vectors for F^. 

Remarks. 
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(1) Note the drop in regularity: we require C^^" initial data for C^^" (in 
space) solutions. This occurs also in the main result in [BaconneauLunardi] . 
and for the same reason. It was also observed in [Freirej . where local exis- 
tence for the 'one-phase' problem is obtained by a different argument. 

(2) The result, at this point, lacks both a uniqueness statement and a 
continuation criterion. It is conceivable that assuming higher regularity on 
the initial data (C^"'""), a uniqueness statement could be obtained following 
the argument in |BaconneauLunardi| . Both issues will be considered in 
further work. 

(3) The restriction to two dimensions makes the argument technically 
easier, but is probably inessential. On the other hand, generalizing to con- 
figurations beyond graphs may require a completely different line of argu- 
ment. 

Outline. Section 1: We derive the compatibility condition and the equa- 
tion of motion for the junction A(t) (and for the 'interface' Tt). Section 2: 
In a standard way, we use a time-dependent diffeomorphism to pass to a 
system (with four unknown functions) over the fixed domains Dq-. 

w\y, t), y e dI — > u\x, t),p{x, t), x E Dq, 

where p is an extension to of the function pt parametrizing Tt as a normal 
graph over Tq. Section 3: Following the argument in [BaconneauLunardij . 
we introduce a transformation of the dependent variables: 

u\x, t),p{x, t) — > U\x, t), X £ D^, 

which leads to a system on x [0, T] with three unknown functions and three 
conjugation conditions on Fq. Section 4: We describe the fixed-point scheme, 
the associated linear problem and the required Holder space estimates. This 
concludes the proof, except for verifying that the linear system satisfies the 
'complementarity conditions', which is done in Section 5. 
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at a seminar at the FU Berlin (June 2007), and the current version at the 
conference Nonlinear PDEs at IMPA (August 2008) . Finally, it is a pleasure 
to thank N.Alikakos (University of Athens) for originally suggesting the 
problem and for his interest in the results. 
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1. Compatibility conditions. 

We wish to find solutions that converge to the initial data in a strong 
sense, as f ^ 0+. It is well-known from parabolic theory that, for such 
solutions to exist, the initial data have to satisfy compatibility conditions; 
in the present case, of 'order 0' and 'order 1'. Before going any further, 
one should check that the geometric problem is compatible with a graph 
formulation, in the sense that one can find reasonably large families of initial 
configurations satisfying these compatibility requirements. In this section 
we state the conditions, analyze their geometry and examine their internal 
consistency. 

'Order 0' compatibility is just the requirement that the initial data sat- 
isfy the conjugation and boundary conditions (i.e., that the define a 
configuration). 

On Tq this corresponds to two zero-order conditions on the (pointwise 
matching) and two first-order conditions (angle). Assuming the ^t'oiTo 
given (and hence the curve Aq), the conditions on the first-order data (the 
tangent planes to S^(0) along A(0)) are clearly compatible. (Construction: 
at each point, on the plane orthogonal to the tangent line of Aq- which is 
determined by zero-order data- we can arbitrarily pick one unit vector, say 
A^^(O), and then the other two are uniquely determined; so we have one real 
degree of freedom.) 

'Order one' compatibility conditions along Fq are obtained by differenti- 
ating the incidence conditions in t and using the equation, bearing in mind 
that r is time-dependent. 

Parametrizing Tf by a{s, t) € (,s € M/27rZ is a fixed periodic variable, 
not arc length), we compute (using the equation of motion) the velocity 
vector at points of A^, = ■§j:{G^ {a{s,t),t)): 

= [a, Dw^ ■ a + dtw^] = [d, Dw^ ■ a + H% 

where a = dfa. 

The compatibility condition is then the statement that, at t 
fortiori, for all t): 

V^ = V^ = y e Tt. 

In principle, this corresponds to up to six independent scalar conditions. To 
analyze the situation, observe that at points of A^ we have three natural 
orthonormal frames {E, ,T^}, where E is the unit tangent vector to At 
(orienting counterclockwise) and we define = EAN^ (vector product). 



= (and, a 
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a tangent vector to S-^; in particular, T^ + T^ = . Noting that {V^ , N^) = 
we may write: 

Six independent scalar conditions equivalent to equality of the are then 
obtained as follows: 

(1.) Add {V\N^) = {V^N^) and {V^N^) = {V\N^). 

(2.) Add {V\T^) = {V^T^) and {V^,T^) = {V^T^). 

(3.) Take the difference {V^,N^) = {V^, N^) minus {V'^, iV^) = (y3^ ^2jj_ 

(4.) Take the difference {V\T^) = {V^,T^) minus (F^T^) = {V^,T^). 

(5./6.) {V\E) = {V\E) = {V\E). 

As we proceed to write these relations in a more explicit form, we will 
preserve their labels (1)...(6). We remark that conditions (1)...(4) already 
occur for curves, while (5) and (6) are new for surfaces. 

In terms of the coefficients and A''^, we easily obtain (using N^ + N'^ = 

and T^ + T^ = T^): 

{1)H^ + ij2 = H^. 

(Here we used (T^, - N^) = ^/s and {N^,N^ - N^) = 0.) 
(4)/xi - = -^H\ 
(5/6) Ai = A3 and A^ = A^. 

To make these relations explicit, one has to compute the , X^ and 
fi^ , which is elementary. Denote by r the unit tangent vector to Fj, and set 
Ve = \/l + {drW^Y (independent of I on Tt) and (Dw^)-^ = {—d2W^, diw^). 
The results are: 

= ^[n~ {drW^)iDw^)^,dnW^], 

V Ve 

/ = - — ((a • + {dnW^)H^), 
Ve 

X^ = — (d • r + {drw^)id ■ Dw^ + v^H^)). 

With this we can rewrite (l)-(4) as: 
{1)H^ +H'^ = H^. 

(2) (d • n){v^ + + {dnW^)H^ + {dnw'^)H^ = (d • n)v^ + {dnW^)H^. 

(3) H'-H^ = {V^/vE){{d ■ n)v' + {dnW^)H^). 

(4) (d • n){v^ - v'^) + {dnW^)H^ - {dnW^)H^ = -VSveH^. 
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When rewriting conditions (5/6), we observe that terms involving the 
tangential component a • r of the moving boundary velocity cancel, and 
there is a common factor drW = d-rW^ (independent of /). Wc have: 

(5/6) {drw){{a ■ n)dnw^ + v^H^) = {drw){{d ■ n)dnw'^ + v'^H^) 
= {drw){{d ■ n)dnW^ + v^H^). 

Thus (5) and (6) occur only at points where dxw ^ (that is, where is 
not 'horizontal'). Regarding the zero and first order data of the configuration 
as given, (l)-(6) give six homogeneous linear relations for the four 'second- 
order quantities' a-n and , so it is conceivable that only the zero solution 
might exist. That this is not so is part of the conclusion of the following 
proposition. 

Proposition 1.1. The homogeneous system of linear relations (l)-(6) 
has rank two, and is generated by (1) and (3). 

This has the following interpretation: given zero and first-order data, 
the only constraint on the second-order data of an initial configuration of 
surfaces satisfying 'order one compatibility at t = 0' is condition (1) on the 
mean curvatures at the junction curve, while (3) specifies how the mean 
curvatures determine the normal velocity of the interface Ft, in particular 
at t = (since the normal velocity d-n is independent of how we parametrize 
Ft, we denote it by r„): 

Geometrically in R^, condition (3) says that the component y"""" (normal 
to the junction Aj) of the velocity vector, which is independent of how one 
parametrizes A^, is determined by the via: 

v3 

Proof. This is computational and elementary, so we just describe the 
main steps. One just has to observe that (as for curves), the constraints on 
first derivatives arising from writing + iV^ = componentwise can be 
'solved', in the sense of the following elementary algebraic lemma: 
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Lemma 1.2. Given real numbers a > and a^,a^, a^, let = 
suppose these numbers satisfy the following relations: 

1 1 1 




V V V' V V v 




Then (assuming > to normalize the labeling): 




Remark. This lemma has the following equivalent geometric formulation: 
let u>^,uj'^,uj^ be vectors in M? of the same length, satisfying a;-^ + = ui^. 
Then, denoting by Ro the counterclockwise rotation operator in (by 6 
radians), we have: 



Since n is the inner unit normal of Dl = D|, the condition > 
corresponds to labeling the surfaces so that is below T,^ near A^. 

To prove proposition 1.1, we substitute the values of and dnW^ /v^ 
given by lemma 1.2 (for / = 1,2, in terms of the same quantities for / = 3) 
into relations (2)-(6). Using relation (1), we find that each of (2) and (4) 
is equivalent to (3). Where drW / 0, adding and subtracting the relations 
obtained from (5) and (6) with this substitution, we recover respectively (1) 
and (3), completing the proof. 

There is no 'order 1' compatibility condition to consider at the fixed 
boundary dCl, since the Neumann condition d^w^ = is first-order. 

2. An equivalent system on a fixed domain. 

In this section we describe new dependent variables {x,t), p{x,t) de- 
fined on fixed cylindrical domains Dq x [0, T], J7 x [0, T] (resp.), and a 
parabolic system with conjugation conditions, equivalent to the original sys- 
tem for w^y,t) on G [0,r]. 

2.1. An extension operator. Fq C $7 C is a simple closed curve 
of class C^+", bounding domains Dq = Dq on the inside and Dq = Q — Dq; 



a;2 = R5^/s[u% 
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denote by v the inner unit normal vector field of Fq (of class C^"*"") and by 
r : A/q M the oriented distance to Fq, where A/q is a fixed tubular neigh- 
borhood of Fq. We have r G C'^^"(7Vo), and the natural extension v = Dr 
of ly to A/q- With f := sup^v/^ r, consider the smaller tubular neighborhood: 

X:={xe A/'o; |r(a;)| < ro := min{f/2, 1/2}}. 

Denoting by ^ S C°°(M; [0, 1]) a fixed function supported in [— ro,ro], equal 
to 1 in [— ro/2,ro/2] with \('\ < 3/ro everywhere, we see that C{^)i defined 
as (^or in AA and identically zero outside of AA, is in C^^°'{M?) and coincides 
with r in M := r"^(-ro/2, ro/2) c J\f. 

Using we extend v = Dr from A/i to M?, by setting := ^z/ G 
C2+°(M2;M2) i7 is supported in M and satisfies |p| < 1 in R^. Note that 
D,yU = in A/i, while: 

which is supported in J\f — A/i . 

We will need a fixed extension operator C^^°'{Tq) C'^'+"(]R^), defined 
in a standard way using C (for k>2). Given p G C''+"(Fo), extend p to A/'o 
by projection along normal line segments (so di,p = in A/'o), then set: 

r Cp, Mo 
^' \ 0, xeR^-M 

It is clear that p is supported in AA, and that ||p||c*=+«(M2) ^ c^l|p|lc'=+«(ro)' 
where c^; depends only on C and on | |r| |(^3+a(jv-g) . Since the extension operator 
is ^-independent, the same estimate holds for parabolic Holder norms when 
p depends on t. Note also: 

Dp{pv) = C,v{p)y + pDpi? 

= ({{d,0p + 0)u + p{d.0iy 
= 2{d,0pi>, 

supported in AA — A/i . 

If U is already defined on a set containing Fo, we denote by £[U] the 
restriction-extension operator: 
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2.2. Diffeomorphism. For t > small enough, the interface Tt will 
be contained in A/i, and we may parametrize it as a normal graph over Fq: 

Tt = {x + p{x, t)i^{x);x G To}, for some p : Fq x [0, T] R. 

Recall the subscript t always denotes 'function (or map, or set) at time t\ so 
Ptix) := p(x, t). We use the extension operator £ to define a diffeomorphism 
ift of R2 (or of n), for t G [0, T]: 

(pt{x) = X + pt{x)u{x), (po = id. 

For t e [0, T] small enough, ipt is a diffeomorphism (of class C^+"^), equal to 
the identity in M? — M, and mapping: 

n^n, Fo^Ft, Di^Dl 7 = 1,2,3, 

diffeomorphically in each case. If p G C^+"^'^+"/^(Fo x [0, T]), is in the 
same Holder class as a function of {x,t). The differential of (pt is given by: 

D^t = I + {Dpt)i) + ptDy\ 

in particular, in light of the calculation in 2.1, Dipt maps P via: 

£>¥)t[p] =v + Di,{ptu) = [1 + 2(d^C)ft]i^, 

which implies: 

Here C(p) e C^+"(M?), and is identically 1 in A/l n (R^ _ AT), with modulus 
bounded above by 1 everywhere (provided only we have \pt\ < ro/12 for 
t G [0, T]). In particular, on Fq: 

(recall u := Dr on A/o). Of course, z/(x + pt(a;)i/(x)) is in general not normal 
to Ft. 

2.3 Equations over Dq. Assume : ^ R are solutions of 
(MCM)^ for t G [0,T]. Define, for I = 1,2,3: 

u^{x,t) =w^{ipt{x),t),x G D^; w\y,t)=u\il;t{y),t),y G 
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(where il>t '■= {'Pt)~^)- A standard calculation yields (omitting the super- 
script / sometimes, and using the Einstein summation convention through- 
out): 

= dtw- g'\Dw)Dl^w = dtu - h'^'^Dl^u + idcu){dtr - g'^Df^r), 

or equivalently: 

dtu - h'^'Dl^u - {d.uMrM'p' - h'^'dl,^^) = 0, 

where hab = da^^d},<^gij is the pullback metric (of g under </?). Since the 
parametrization over Dq of is: 

Fl{x) = Gl{ipt{x)) = [pt{x),wl{^t{x))] = [^t{x),ul{x)] G X R, 

we have for x G D^: 

hi, = {daFl, dbFl) = da<pt ■ dm + (dauOidbui) 
= Sab + da{pt^) ■ dbiptv) + {daul){dbul) 

■■= Sab + KbiPt, Dpt, Dul). 
Define the operator on functions (or maps) f\x,t) on Dq x [0, T]: 

The equation on Dq is: 

We may express L^lf] in terms of p and u: 

Lh[p\ = - ptrhD^u - 2{Dp, Du)h. 

Using {Dlp)^'^[u] = C,{p)^: we find: 

{D^)-\Lh^] = ap)LhW - {D^)-\ptrhD''v + 2{Dp,Du)h]. 
Hence the system in the variables {u^,p) is: 

V[«VC(p)(dptx0i>/,/[p]+W-(I)^t)"i[ptr;,.I)2i^+2(I)p,I)z7)^,] = 0, (PDi?)(. 
where: 

hib = Sab + h'abip, Dp, Du^), Dipt=I + D{pu), 
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{Dp,Di>)f,i = ih'r\dap)id,i>)e 



2.4. Boundary conditions at Fq. (More precisely, these are 'conju- 
gation conditions', since Fq is an 'internal boundary' in $7). We begin by 
computing how Dipt^Y'o ^'^^^ orthonormal frame {i^,to}, tq := —ly^ 

(where _L denotes counterclockwise 7r/2 rotation in M^). Parametrizing Fq 
by arc length s, we have: 

ro(s) = To, ly' = -koTo 

(with ko the curvature of Fq). Composing with (ft, we obtain the parametriza- 
tion of F^: Ft(s) = Fo(s)-)-pt(Fo(s))^'(s), with tangent and inner unit normal 
vectors: 



F((s) = (1 - kop)To + idroP)i^, n 



(1 - kop)l^ - (rfroP)To 

[{i-kopr+dr,pn/^' 



Using Di,u = 0,di,p = 0,DtqU = u' on Fq, we find for an arbitrary v G 
and X G Fq: 

Dipt{x)[v] =V + {V ■ To)idrop)l' - (V ■ To){kop)To, 

in particular verifying again: 

Diftlv] = V, DiptiM = (1 - kQp)TQ + {droP)v = Fj(s). 

We need vectors mapping to n and to r := —n^ under Dipt. It is easy to 
see that, defining: 

To := [(1 - k^pf + {dr^,pf]~'^^'^TQ != anTQ, 

/. = (!- M"H[(1 - M' + {dr,pf]^'^u - «p)fo} 

we have for these vector fields p, = p{p, dj-gp), tq = To(p, dr^p)'- 

D(pt{x)[p] = n, Dipt{x)[fo] =T, a; e Fq. 
(The expressions with the are used in section 4.) Hence, for x G Fq: 
dnw\ipt{x)) = d^j,u\x); drw\ipt{x)) = dfou\x), 



13 



V' = [1 + {drw'f + {dnw'ff^ = [1 + {dr.u'f + [d^f]^'^ 

:= g{p, droP, drgU^, dyU^). 

{Remark: It is easy to see that /x — z/ = — ((iT-o/o)ro + 0(p^ + IdroPp); at t = 0, 
11 = V on 7o, since p = 0.) 

We now state the conjugation conditions in terms of ,p'. 

u\x, t) = u'ix, t) = n3(x, t) {BC'^)(u',p) 
^ + ^ = ^ {BC^)(^ui,p) 

In the computation of the second angle condition, we use the fact that dr^ 
is independent of / on Tq (which follows from the matching condition). Here 
= Q{p, drgp, d^-gU^ , di,u^) has the expression given above. 

Since the diffeomorphism ipt is the identity on dQ, we have the additional 
(Neumann) boundary condition for u'^: 

d^u^ = on dn. 

(Recah dD^ = Tq U Sri.) 

3. The transformed system. 

The system described in the preceding section includes 4 scalar conju- 
gation conditions, but only three evolution equations for the four unknowns 
, p; there is no explicit evolution equation for p (or its extension p). Yet 
(assuming the evolution of triple junctions is well-posed as a PDE) it is clear 
geometrically that the evolution of the determines that of p. In this sec- 
tion, following the technique introduced in the paper |BaconneauLunardi| 
we show it is possible to define new dependent variables on x [0, T] 
in such a way that p is recovered from the via the extension operator £. 
The drawback (as in [BaconneauLunardij ) is having to introduce 'non-local 
highest-order terms' in the resulting equation for the . 

Following |BaconneauLunar"di] . define : Dq x [0,T] ^ M via: 
u^{x, t) = Uq{x) + {duul){x)p{x, t) + U\x, t). 
(Note C^i^^Q = 0.) The matching conditions {BC^)(^^i on Fq imply: 

pdyul + = pduul + = pdyul + on Fq. 
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Equivalently, we have the following equalities on Tq: 



P\ro 



difUQ di/ZiQ di/VjQ di/VjQ di/U^ di/ZiQ 
The last equality may be written in the form: 

— U'^ dyUQ — dyUQ 

— dyll^ — dyUQ Vq ' 

which we regard as a new form of the matching condition. We used the fact 
that the angle conditions at t = imply: 

\{dT,u\ — dj,ul) = \{duul - duul). 

^0 ^0 
This new matching condition may also be written in the form: 

^3 _jj2 _ jjl 

2 + 1 = 0' 

or equivalently (using the first angle condition at t = 0) : 

f/l 7t2 7t3 

— + — = — onTo {BC%i, 
K ^0 ^0 

which we adopt as the matching condition for the system. 
The main point of the method is that the relation: 

p = 5o{U^ - C/i) on To, := {duu}, - d^ul)-^ 

allows us to express p in terms of £\U'^ — U^\- 

To explain how this is done, we take a moment to examine the extension 
operator applied to the product of two functions, /, g defined on Fq. Denote 
by the subscript rad the extension to A/q constant along normals of a function 
defined on Vq. By definition, on A/q we have: 

Ua) = Cifg)rad = if)rad9 

(this is true on A^, and both sides vanish on A/q — AA). Let x be a smooth 
function in M^, equal to 1 on A^ and vanishing outside A/q- Given / on Tq, 
define / on by: 

/ = Xfrad on Mo, / = on - Mo- 
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Then it is easy to check that: 

{fgj=fg on M^. 



Thus we may recover p from — as fohows: 

p = 6o£[U^ -U']. 

We now express the angle boundary conditions in terms of the . Since: 
= ui + SoiW^ - U^)d^ul + on Tq, 

we have on Tq: 

droU^ = droUQ + droiSoiU'^ - U^)duUQ] + droU^ 
:= AtaniU^ - U\ dr,{U^ - U^)) + d^^U' 

{Atan for 'afHne', as a function of the indicated arguments.) For the normal 
derivatives, we use the expression defining on Dq, and recall that d^p = 
on Tq; hence (using also Di,v = Qi): 

duu\p^ = d^ui + (5o([/2 - U^)D\i{iy, u) + d^^ 

:= Ar^oriU^ - U^) + d,U' . 

In section 2 we found functional expressions = G{p, dr^p, d^-^^u^ , dyU^). If 
in the expressions defining the vector fields tq and p on Fq we replace p and 
droP by 5o(f^^ — U^) and drQ[5Q{U'^ — U^)] (resp.) and use the expressions 
just found for dT^u^ and d^u^ , we obtain: 

vl=QmU^-U^),dr,[6o{U'^-U%Atan{U^-U\dr,{U'-U'))^^^ 



From now on we adopt the notational convention for the summation of 
quantities depending on a superscript / = 1, 2, 3: 



:= o} + - 
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It follows from the above (using also the boundary conditions for the Uq) 
that the boundary conditions take the form: 



I e; ^ + e; ^ + so{u' - u') e; = {Bc'h' 

Computing the equation satisfied by the is straightforward. First we 
write down the expression for Lf^i[u^], isolating the lower-order terms: 

where: 

6l}\~p,D~p) ■= -2{Dp,D{dWo))^i -trf^iD^Q - ptr^iD\d^ui). 
Now using the equation [PDE)^^i p^ from section 2 and the definitions: 
Cf,iip,Dp) := {D<p)-^[ptr,,iD^u + 2{Dp,Di?),,r] with = I + L>(pP), 

Dp) := D{ui + pdpui) ■ Chi, 

we find: 

L^ip'] + [{d^ui - + pd^ui + U')\Lf,,[p\ 

+DU' ■ ey{~p, Dp) + 6^}\p, Dp) + Cf]\~p, Dp) = 0. 

(Note that {Dip)~^ depends on Dp, and that the first-order terms involving 
the metric depend on DU^ , in addition to p and Dp.) 

This may be regarded as a quasilinear system in C/^, if we agree to replace 
p by 5q£\U'^ — U'^] at every occurrence of p or Dp. Due to the presence of 
the restriction-extension operator this introduces 'non-local terms', even 
to highest order. 

To make the dependence on a bit more explicit, compute: 

= SoLhi[£[U^ - U']] - 2{D5o,D{S[U^ - U']))^! - S[U^ - U^]trf,iDH^o. 
Thus we have the equation for U^: 

Lh,[U']+A{p,Dp,DU\^,~^,^^,_^,^L^,,[£[U^ - U^]] 
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where: 

A{p, Dp, DU^) = 5o[{dpui - ap)dv{ui + pdpui + U^)]; 

T(j = 2{D6o,D{£[U^ - U^]))f^i +£[U^ - U^]tr^^iD^6o 

-D{ul + pdpM^ + U^) ■ Chi 

+2{Dp, D{dpul))hi + tr^iD^ui + ptr^i D^{dpul), 

and Cf^i was given above. 

We record this in abbreviated form (changing the notation sUghtly to 
exhibit the dependence of and ^ on U): 

The critical term is AIjL^i^ [<5[{7^— f/^]], due to the 'non-local dependence' 

on (i-e., it is not given by the action of a differential operator on the , 
so we may not quote standard linear parabolic theory at this point). This 
term has the important property: 

4=0 = on Dl 

This is easy to see from the definition, since is 5q multiplied by the 
expression: 

djjul - C{p)dp{uQ + pdpul + U^). 

At t = we have = 0, p = (since p = 0) and therefore C(p) = 1- This 
makes it possible to treat this linear system as a small perturbation of a 
standard linear parabolic system, as long as T is small. This was also the 
general strategy adopted in [BaconneauLunar"dI] . but we actually deal with 
this term in a different way (in a sense, less sophisticated; see Section 4). 
P\t=o = also implies C^i = at t = 0. 

Finally, we record the dependence of the metric on DU^: 
hib = Sab + daipi^) • dbipi^) 

+ {daui + daipdpui) + daU'Kdbui + dbipdpui) + dbU'), 

where we set p = do£[U'^ ~ U^]- 
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Note also the two properties of the term JT^: it involves derivatives of 
up to first order, and equals tVf^iD'^UQ = VqHq at f = 0. 

4. The associated linear system. In this section we describe the 
(standard) fixed-point argument used to prove local existence, based on the 
solution of an associated linear system with 'non-local terms'. 

4.1 Linearization of the conjugation conditions. 

The angle conjugation conditions {BC)ui involve sums of so we 

begin by considering the linearization C oi 1/v sA, points of Fq, where: 



■"iro = V 1 + {d^,uY + {dfouy = v{p, droP, DU) 
and we linearize at p = 0,U = 0. We find: 

£{-) = — \[{duUo)C{dfj,u) + {droUo)C{dfou)] 

V Vq 

Prom the expressions for tq and /x given in section 2: 

To = aiiTo, P = a2lTQ + a22V, (^ij \p=Q = ^ij ^ 

we have: 

C{dfou) = droU + {di,uo)droP + C{au)droUQ + d^uo(ro, i')p, 

C{d^u) = dvU + C{a2i)dr„UQ + C{a22)dyUQ + d^uo{u, v)p. 

Straightforward computation from the expressions for the a^j given in section 
2 shows that: 

L{aii) = kop, C{a2i) = -droP, ^{022) = 0. 

Substituting in the expression for C{v~^), we find that the coefficient of dj-oP 
vanishes, yielding: 

£(-) = — \{Duo ■ DU + [di.uod'^uoii', v) + droUod'^UQ{u,To) + {droUQ)^kQ]p} 

V Vq 

Restoring the superscript /, we write this (always at Fq) in the form: 

^ = ^ - TT^^^o • DU' + 7oV + e(-'^\p, droP, DU') 
'"0 
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where: 

^0 = -7--]T^[duui(ful{i^, v) + droui(ful{iy,To) + [drouifko] 

K) 

and = 0{p^ + {dT^pY + \DU^\^), with constants depending only on 

the initial data. 

The first angle condition {BC^)ui then becomes, with E'-^'> = Y^'jE^^^^: 

Now make the substitution p — (5o(C/^ — ^7^), valid on Fq. We find: 

. DU' + jj,'\u' - U') = E('Hp,dr,p,DU)^^=So(u-_u^), 

where: 



The second angle condition {BC^)ui has the form: 

I ^ I ^ I ^0 

or (using the linearization of computed above): 

where \E^'^\p,drop., DU)\ = 0{p'^ + {dr^pf + \DU\^). Again making the 
substitution p — 5q{U^ — U^), we find: 

^'{^d.U'~^Dui-DU')+ji,'\u'^U') = E('\p,dr,p,DU)\^=Soim_ui^, 



I "0 



where: 
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In summary, the conjugation conditions may be written in linearized 
form: 



(with zero initial conditions.) Here we defined: 



^0 



Observe (for future use) that, since 1 — (di/tio)^/(vo)^ = [1 + (^^to^o)^]/(^o)^) 
we have: 



^ . DU' = [1 + {d^,u,f\Y;^-\-d^U' - (d.o^o)E'^^^o^'- 

(The notation uses the fact that (iro^o independent of / at points of Fq). 
4.2 The fixed-point scheme. 

We write the quasilinear system {PDE)iji / [BC)ui in a shghtly modified 
'hnearized form', as follows: 

{LPDEILBOv I ^, ^ ^ jg^^^^j ^ ^^^^j ^ j.^^ ^^^3^^ ^ ^ 

where: 

/^(x,t) := .F/, + [(/i^)^^ - (/.f,)'^'']<,c/^ - - v^W 

and we collapsed two linearized angle conditions into one, introducing a 
slight change in notation. h\ = is the induced metric at t = 0. 

This will be solved by a standard fixed-point argument in the space: 

Xi?,T := {U = {U\ U\ U'); U' G C^+a.iWs^^/ ^ [q, T]), ^,{^0 ^ 0}, 

with suitable choices of R and T, to be described soon. (We organize the 
fixed-point argument in the same way as |BaconneauLunardi| . ) 
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Define an operator ^ : V U hy assigning to F G Hr^t the solution of 
the linear system: 

,^p^p J L,,^[U^]=fU^,t) onDlx{0,T] 

^ °^ \ E/^ = 0, Mo[DU] + 7o[C/] = Ev on Tq, d^Uf^^ = 

Consider the following assumptions for system (LPDEq): 

(1) (Regularity of coefficients) G C"(£»^), Bo,7o G C^+"(ro), (i;^)"^ G 
C2+"(ro); 

(2) (Complementarity) The conjugation operator {Yli'^:^o[DU]) sat- 
isfies the 'complementarity conditions' on Tq, with respect to the operators 

Under these assumptions, it is a classical fact that this system (with 
vanishing initial data) has a unique solution, provided fy G C"'"/^(Dq x 
[0,r] and Ev G C^+^'^^+^^^iTo x [0,r]) satisfy also: 

(3) (Compatibility with f/|^^o = 0:) 



0. 



/ '"o |t=0 



The solution satisfies the estimate: 

||C/'||2+a<Mo(^||/f|U + ||Ey|| 
I 

with Mq depending only on the Holder norms of the coefficients in (1). 

The verification of condition (1) under the assumption Uq G C^~^'^{Dq) 
is straightforward. Verifying (3) is also easy, since when = wc have 
Eu = and = tr^iD^UQ, so that the compatibility condition amounts 
to: 

Y^^tr^iDW, = 0, 

which is just the compatibility condition considered in section 1 (equivalent 

to E/^o = on To.) 

Complementarity (assumption (2)) also holds for the system (LPDEq), 
but verifying this is more technical- it is done in the next section. Assuming 
(2) for the moment, we conclude the map $ is well-defined. 

4.3 Contraction estimates in Holder norms. 
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To finish the argument, we must verify that (with suitable choices of R 
and r) $ maps into Xt,r and is a contraction in this space. 

Let Vi,V2e Xr^t- Then W = ^{Vi) - ^{V2) is the solution of the linear 
problem: 

(TPDF J ^'^^[^'1=^^,-/4 onDo^x(0,T] 

^ \ E'/^ = 0, Mo[DW] + jo[W] = Ev, - Ev, on To, d^Wf^^ = 

(with zero initial conditions.) Thus we have the estimate: 

||W^||2+a < Mo{J2\\fv, - /4lU + ll^^^i - ^y.lll+a). 

We write the difference fy^ — fy^ as the sum of five terms: 

j^i I'l 'T7J 

JVi ^ JV2 — -^Vi ^ •^V2 

+{hit - h\;^)dl,Vl + (/i^ - h\;^)dl,{Vl - Vi) 
+(^(.i-OV^Jf(V/)] 

+^(..(/ii.f -4f)9>/ 
+</ii.f 9i(y/ - Vi) 

.- f{^)I + . . . + . 

We need bounds for each of the Z*-*^^ in terms of 1 1 Vi — V2I |2+a, R and T. 
The details are lengthy but standard, and it suffices to state the estimates 
with a brief justification, (co denotes a constant depending only on the initial 
data, which may change from one occurrence to the next.) 



11^4 - •^4lU < Cor°/2||Fi - V2\\2+a, 

since J^y depends only on V,DV,S{V) and D{6{V)), and vanishes at t = 
0. (The 'nonlocality' of £ docs not make the estimates harder, since this 
operator is linear in V and bounded in (72+a,i+a/2 j^orm.) 

< CoT-^^m - V2\\2+aR + CoT"^^\\Vi - V2\\2+a, 

since hy^ — hy^and — hy^ depend only on derivatives of Vi, V2 up to first 
order and vanish at t = 0. 
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(^(.1 - OL^^^ [£{V-[)\ < CoT-/^\\V, - V2\\2+aR, 

since Ay^ — Ay^ depends on derivatives of Vi , V2 only up to first order and 
vanishes at t = 0. 

since Ay^ depends on derivatives of V2 up to first order, and vanishes at 
t = 0. 

Finally, 

AM'dlbiVi' - Vi) < coT^/^IIFi - V2\WaR\ 

for the same reason. 

A bit more involved (but also needed) is be the estimate over Tq: 

\\Ey, - £;yj|l+a < CoT^/'lll^l - V2\\2+aR. 

Although the 1 + a norm of Ey involves D^F, the estimate holds since Ey 
is of quadratic order in V and DV . 

We conclude that $ is a 1/2 contraction, provided we choose R and T 
so that 16co(l + i?2)r"/2 < i (16 = 5 X 3 + 1). Assuming this inequality 
holds, applying the contraction estimate to the case V2 = yields, for any 

V G Xr,t: 

||$(V^)||2+a<||f(0)||2+a + ^i?, 

where $(0) is the solution (with zero initial data) of the system (LPDE)q 
with data fy = try^i^D^UQ, Ey = 0, satisfying the estimate: 

mO)\Wa<Mo\\tr^,D''ui\\a. 



The choices of R and T are made as follows. As observed earlier, if the 
geometry of the configuration at time t is sufficiently close to that at 
i = 0, all the geometric constructions used make sense; in particular the 
operator £ and arc well defined. Choose r > so that 1 1 1 1 (72 < r 
quantifies sufficiently close to Uq' . Then pick Tq > so that C/|t=o = 
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and ||C/||2+a < To imply | |(72(£)i) < r. Now choose R > sufficiently 
large to ensure: 

MoWtTf^iyuiWa, < ^R. 

Finally, given R we pick T < Tq small enoug h that 16co(l + i?2)r"/2 < i 
holds. 

With these choices we guarantee, on the one hand, that <^ is a contrac- 
tion; and, on the other, that: 

\mV)\\2+a<\\m\\2+a + lR<R 

for any V G ^r,t- Thus $ maps into ^r^t, and has a unique fixed point 
U. This fixed point is the unique solution of {LPDE)i/ /{LBC)u with 
zero initial data. This concludes the proof of local existence for the sys- 
tem {PDE)^^i p-^, and hence for the original system {PDE)^i. (Except for 
the verification of complementarity for {LPDE)q, carried out in the next 
section). 



5. The complementarity condition. 

In this section we verify that the system of conjugation conditions along 
Fq satisfies the 'complementarity conditions' (Lopatinski-Shapiro) with re- 
spect to the 3x3 linear parabolic operator Lgi^[U^]. (We use the conditions 
as stated in jEidelmanZhitarasuj .) Explicitly, the operator is: 

LgiU' = dtU' - gi'^'dl.U' on 17 X [0,r], 



Sab + daUQdbUQ, 



lab 

9o 



Jab 



dau^dbul 



K) 



Vq = \Jl + (drouly + {dyU^Y on Fq 
The conjugation operator has three components: 



Bi[Y] = llj^.d^U' + XoEi'j^dr,U\ 
^ B2[Y] = (1 + Xl)J2/^dM' - AoE/'(^<C^', 

where in this section we adopt the notation: 



Ao :— dro'Woiro- 
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Fix xo G To and 'straighten the boundary' via a difFeomorphism x from a 
neighborhood of G (with coordinates z = {zi, Z2) to a neighborhood of 
Xo in ri, and with the mapping properties: 

X : {zi = 0} ^ To, x(0) = Xo, dx(0) : d^-, t-^ i^{xo), dx{0) : d^^ h^- to{xo). 

Let Y^{z,t) = {x{z),t). From this point on the symbols 7^, Uq, Aq, 
Vq will denote the values of the corresponding functions at the fixed point 
Xo G To. 

The transformed operator is: 

Lo[Y'] ■.= dtY'-j''^dl,.Y'. 
1^ = x(0)*fi'o(^o) is the pullback metric tensor, with components: 
7ii = 1 + {nif^ 7i2 = -^0^0 , 722 = 1 + -^o> 

and inverse: 



1 + -^0 112 _ Aon^ ,,722 _ 1 + (^0)^ 



0^ 



In the new coordinates, the components of the conjugation operator are: 

Bi[Y\ = Z/^d,,Yi + XoE/j^d^^y, 
^ B2[Y] = (1 + Xl)J:i'^d,,Y' - \,Y.i^.dz,Y', 

The next step is to consider the Fourier-Laplace transform y^(r, 
of Y^ {z\, Z2,t): Fourier transform in Z2 corresponds to the variable ^ G M, 
while Laplace transform in t corresponds to p G C. Furthermore, we adjust 
the signs so that {r > 0} corresponds to = Dq (that is, t = zi for Y^ 
and Y^) and also to -Dg (that is, r = —zi for 1"^). (This is how we deal with 
the fact that we have a conjugation problem, rather than a boundary-value 
problem.) This introduces a sign in the transformed operator 
1,(7^ = —1), which is an ordinary differential operator in the variable r, for 
each fixed {p,C)'- 



dr"^ dr 
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We are interested in solutions of the corresponding ODE which decay as 
r — > +00. If we assume the exponential form: 

this corresponds to lm{pi) > 0. The p/ G C depend on ^ and p, and are 
roots of the indicial equation: 



Using the expressions given earlier for the 7^*-', we find for the roots: 
1 



PI = -nii-^'^'''^ + det[7^]-i + pj^^^ 



1 i- Vq 



where A = + p(l + Ag) and we take the branch Re{y/A) > 0. The 
boundary operators for the ODE are obtained from the correspondence: 

ffyl 

We find: 

Bo[Y] E/^ 

Bi[Y] j:/[^a\ipj) + ^(ie)]y^(o) 

The complementarity condition is the statement that, for each {p,$,) in 
the set: 

A := {ip,0 G C X M; IpI + |^| > 0, (1 + A§)Me(p) > -C'} 

the map ^"^(0) 1-^ defines an isomorphism of C^. In matrix form (after 

canceling a common factor (fg)^^ in each column) we are interested in: 



1 1 -1 
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Lemma 5.1. For each {p,^) G A, det(T(p^)) ^ 0. 

Proof. We begin by introducing some notation. Consider the vectors 
ti G R^,ri e C^: 

Then: 

ripPi ^ Xq ^ _ 1 1 + Aq pi nlXoC _ ± 

~T~ r~r ~ / , ■ .) '-1 ■ '1) i 1 ~ r — "-i ■ 'i > 

and analogously for the other two columns. (Here we use the 'real dot 
product', i.e. tj ■ rj = tj ■ Re{ri) + itj ■ Im(r/)). 

With this notation, we find (after cancelation of common factors and 
using tj ■ rj- = —tf ■ rj) we must show that the matrix: 

1 1 -1 " 

h -n t2- r2 h ■ ra 
_tj; -ri ■ r2 -rs _ 

has non-zero determinant det. Expanding, we have: 

det = (tfri)(t2T2)-(t^T2)(trri)-(iiTi)(t^T3) + (ifri)(i3T3) + (i2-r-2)(ifr3)-(t^T2)(i3T3). 



The argument has two parts. In the first part, wc use the symmetries in 
the conjugation conditions, encoded in the properties of i/: 

tl+t2=t3, \tl\ = l. 

From the geometric version of lemma 1.2, we can express ti and t2 in terms 
of is: 

h - ^n/sh - + -^^3 ' *2 - ^5n/3h - ^^3 ■ 

Substituting these values in the expression for det, we obtain: 

a/3 1 
det = ^ (n + r2) • [{h ■r3)t3 + (4 • r3)i^] + - (ri - rs) • [{4 ■ r3)t3 - ■ 

+r2 ■ {^[(n • t3)t3 + (n • 4)4] + ^[{4 ■ n)t3 - ■ ri)4]}, 
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which simpUfies to give an expression in terms of the rj only: 

\/3 1 
det = — (ri • rs + r2 • rs + ri • r2) + -^{n ■ r2 + r2 ■ r^- + ra ■ rj*"). 



In the second part of the argument, we use symmetries in the rj. Re- 
calhng the expressions for the roots pi, we have (with aj := Me(r/),6/ := 
lm{ri), both in R^): 

r, = a, + »„ „, = A„{l-^-J,-jl, (,, = -^[1,0]. 

We see that 6/ := 6 is independent of /, while the ai satisfy: 



^07 = 0, \ai\ 



a 



|2 



(independent of /). Consider first the imaginary part of det. We have: 

Im(r/ • rj) = Me(r7)Im(rj) + Me(rj)Im(r/) = (a/ + aj) • 6, 

and from this it follows easily that lm{r\ • r3 + r2 • ra + ri • = 0. Similarly, 
since rj- = aj- + ib^, we have: 

lm{ri ■ rj ) = a/ • 6"*" + aj • 6, 

and again we have: Im(ri ■ r2 + r2 ■ + ■ r^) = 0. Thus Im(det)= 0. 
Turning to the real part, since Im(r/) • Im(r/)-'" = 6 • 6"*" = 0: 

Me(r/ • rj) = aj ■ aj — |6p, Me(r7 • rj ) = a/ • a j — 6 • 6"*" = a/ • Oj. 

Using the version of lemma 1.2 for the case J2i — \'^A — I'*!' 

ai = -R47r/3a3, 02 = R2tt/3C'3 

we see that: 

1, ,2 

ai • as = 02 • as = ai • a2 = - -|a| , 

± ± ± V^i |2 

01-02 =02-03 = O3 - Oi = - 
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We conclude: 



det = %-l\a? - m') + l^-^Wf = -(3V3/2)|6p = -(3^3/2)^. 

This is non-zero for any {p,£,) G A (since Re(A) > in A), as we had to 
show. 

Remark. Note that this calculation depends on the fact that the bound- 
ary conditions (that is, the tj) and the roots of the characteristic equation 
(that is, the r/, and ultimately the coefficients of the operator, which come 
from the induced metric at the junction) satisfy the same type of symmetry. 
That is, complementarity of this particular set of conjugation conditions 
seems to be linked to the fact that we are dealing with mean curvature 
motion; this suggests that local well-posedness may fail for a more general 
parabolic system (with the same conjugation conditions.) 
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